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The mode theory of propagation of electromagnetic waves at extremely low frequen- 
cies (ELF) (1.0 to 3,000 cycles per second) is treated in this paper. Starting with the 
representation of the field as a sum of modes, approximate formulas are presented for the 
attenuation and phase constants. Certain alternate representations of the individual modes 
are mentioned. These are used as a basis for describing the physical behavior of the field 
at large distances from the source, particularly near the antipode of the source. At the 
shorter distances, where the range is comparable to the wavelength, the spherical-earth 
mode series is best transformed to a series involving cylindrical wave functions. This latter 
form is used to evaluate the near field behavior of the various field components. 

The effect of the earth's magnetic field is also evaluated using a quasi-longitudinal 
approximation. In general it is indicated that if the gyrofrequency is less than the effec- 
tive value of the collision frequency, the presence of the earth's magnetic field may be 
neglected for ELF. When this condition is not met the attenuation may be increased 
somewhat. The influence of an inhomogeneous ionosphere is also briefly considered and. 
finally, the propagation of KLF pulses is treated. It is suggested that certain observed 
characteristics of ELF waveforms may be attributed to the inclination of the current 
channel in the lightning discharge. 



1. Introduction 

The propagation of electromagnetic waves at extremely low frequencies (ELF) has 
received considerable attention recently [1 to 5]. 2 Unfortunately, the experimental data in 
this frequency range (1.0 cps to 3,000 cps) is still rather limited [8 to 15]. There is little 
doubt, however, that the bulk of the energy is transferred via a waveguide mode. The 
bounding surfaces of this waveguide are the ground and the lower cd^c of the ionospheric 
E region. In theoretical approaches to this problem, it is usually assumed that the iono- 
sphere is homogeneous in its electrical properties. On comparing theoretical attenuation 
rates based on this assumption, it is apparent that the model must be refined to take account 
of the variation of the electron density with height. Such an extension to the mode theory 
has been made recently [16]. The influence of the earth's magnetic field is also expected to 
be important at ELF since the waves are reflected at a level where the collision frequency 
is of the same order as the gyrofrequency for electrons. The ionosphere is thus expected 
to be anisotropic. Unfortunately, the experimental data is not sufficient to confirm whether 
the anisotropy is significant. 

Another feature of ELF propagation is that the distance between source and observer 
may be comparable to the wavelength. For example, at 300 cps the wavelength is 1,000 
km. Many of the experiments are carried out at distances of this order and the computed 
attenuation rates must be considered with caution. 

It is the purpose of this paper to survey the mode theory as it applies to ELF propaga- 
tion and to develop several extensions. Many of the results are presented in graphical form 
as an aid to application. 



1 Contribution from Central Radio Propagation Laboratory, National Bureau of Standards, Boulder, Colo. Taper presented at 
Conference on the Propagation of ELF Radio Waves, Jan. 25, 1960. 

2 Figures in brackets indicate the literature references at the end of this paper. 
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2. Basic Theoretical Model 

The assumed theoretical model is taken to be a homogeneous conducting spherical earth of 
radius a surrounded by a concentric conducting ionospheric shell of inner radius a-\-h. It is 
convenient to introduce the usual spherical coordinate system (r., 0, </>) . Using this model and 
assuming a vertical electric dipole source located at 0=0, r=a, the expression for the vertical 
electric field, E r , is given by an expression of the following form [4, 16] 

£,=|^SS»^^P,(-co S *) (1) 

2kha z ^o sin vk 

where 

/—average current in the source dipole, 
ds — length of source dipole, 
rj — intrinsic impedance of air space ^1207r ohms, 
P„(— cos 6) — Legendre function of argument —cos 6 and (complex) order v, 

v J r—=ka S n where S n , for n=0, 1, 2, . . ., is determined from the boundary condition 
and is described below, and 

2kh(T n = l l 2 for ™=°> =1 for n=l, 2, 3, 



sin 



2khC n 



The individual terms in the above series correspond to the waveguide modes. In the general 
case, the factor S n is obtained from the solution of a complicated transcendental equation 
which involves spherical Bessel functions of large argument and complex order. Certain as- 
pects of this problem have been discussed by Schumann [1] and also the author [4]. For a 
homogeneous earth and a homogeneous ionosphere, both assumed isotropic, it has been shown 
that S n may be approximated by [17] 



where 



Hs)'i[>+V'+-H]? 



A l _L 1 

A= N+N g 

in terms of the refractive indices, A T g and N i} of the earth and the ionosphere respectively. 
This equation is valid subject to the condition |A|M<<1- The sign of the radical is chosen 
in the above equation which makes the real part positive. The values of S n then are located 
in the fourth quadrant of the complex plane. When n = 0, the above simplifies to 

H^stf" (3) 

Now, since |A|&A<C<0> the radical in eq (1) can be expanded for n^>0 to yield 
Furthermore, if in addition, 

kh ^^ \kh) 
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it is permissible to write, for all n, 

where e =l, e n =2 (for n=l, 2, 3, . . . ). This latter equation is valid when the "waveguide 
modes" are not near cutoff and is in the form given by Schumann [2]. It is quite analogous 
to the standard results for the propagation constant in a rectangular waveguide with finitely 
conducting walls [18, 19]. 

3. Behavior of the Modes for the Spherical Earth 

For purposes of computation of the mode series, several simplifications can be made. The 
asymptotic expansion for the Legendre function, given by 

is valid if M^>>1 and not near or -k [20]. Thus, the modes are simply proportional to 

— -— r cos \kaS n (tt- 0) -j 1 
(si n A)* L 4 J 

which apart from a constant factor can be identified as the linear combination of two peripheral 
waves. These have the form 

0-ikaS n 9 



(sin 0) 
and 

1 



(sin 0)2 
where 0<7r. 

These waves are traveling in opposing directions along the two respective great circle paths 
ad and a(2ir-0) f from the source to the observer. It is noticed that there is a tt/2 phase advance 
which the wave traveling on the long great circle path picks up as it goes through the pole d=w. 
The linear combination of these two traveling waves is to form a standing wave pattern whose 
distance 8 m between minimums is approximately given by 

kb m Re S n =7r or 8 m =\/(2 Re S H ) 

subject to 

Im S„< Re S n . 

The preceding asymptotic expansion for P v ( — cos 0) is not usable at and in the vicinity of 
the pole 6= ir. In this region a suitable representation is given by [20] 

P„(_cos e)=J (v) -f sin s (^) [^- J,(n) +~ </ 3 «]+0 (sin 4 (^)) (7) 

where 77= (2^+ 1) sin ((71— 0)/2). J m (y), for m=0, 1, 2, and 3, is the Bessel function of first type 
of argument 77 and order m. When ir—0 is small, the first term is usually sufficient and further- 
more 
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Thus for mode n, the field in the neighborhood of the pole is proportional to the Bessel function 

Jo[kaS n (w-e)]. 

It is then not surprising to see that the first term of the asymptotic expansion of J is the same 
as that of P v ( — cos 6). 

At the lower end of the ELF band, where the wavelength is large compared to the height 
of the ionospheric reflecting layer, the electric field is essentially radial and only one wave- 
guide-type mode is significant. The field is thus expressed by the first term of the mode series 
which reads 

Er=TTTT2 • -P„(— cos0) (8) 

4:kha 2 sm vir 

where v-^-^^kaSo and 

in terms of the relative refractive indices N t and N g of the homogeneous ionosphere and the 
homogeneous ground, respectively. Now as mentioned above, the factor P v { — cos 6) may be 
replaced by an asymptotic expansion if kad or ka(w—d) is somewhat greater than unity. The 
field in this case may be regarded as two azimuthal-type traveling waves. Furthermore, at the 
pole (0 near tt) where the second of these restrictions is violated, it is possible to use an equivalent 
representation which correctly accounts for the axial focusing. An alternate viewpoint which 
is suitable at ELF is to consider the field as a superposition of cavity-resonator type modes. It 
is expected that such a representation would be very good when the circumference of the earth 
is becoming comparable to the wavelength. A suggestion of this kind was apparently first put 
forth by Schumann. 

The starting point is the expansion formula 

P v (-x) ^ 1 - 2n+l , 

sin™- ~ tt^o nW n(n+l)-v(v+l) K } 

where the summation is over integral values of n. This result follows directly from a formula 
given by Magnus and Oberhettinger [20] (p. 57) which is valid for v^O, ±1, ±2, . . ., and 
O^0<7r. The electric field, for A/a<<X is thus written 



7^+1) £ J*? + * (10) 

where x=cos 6. The early terms of the series are then proportional to 

Po(x) = l, 
Pi 0*0 = cos 6, 
JPi(s) = i(3 cos 2 0-l), 
and so on. 

Retaining just the first term it is seen that 

Er=E r ] n=0 =- — o — r~ (11) 

4ira 2 e h too 

which is independent of 6. Clearly this corresponds to a concentric spherical capacitor ener- 
gized by a current Ids/h resulting in a constant voltage hE° r between the plates. On rewriting 
eq (11) in the form 

hE o = (Ids/h) 

icoC e 
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it is seen that 



C e = 



47ra 2 e 



h 



which can be identified as the capacity between the spherical surfaces whose areas are both 
4wa 2 within the approximation h/a<^<Cl. 

4. Earth-Flattening Approximation 

The possibility that the curvature of the earth may be neglected is now investigated. To 
simplify field calculations at short distances, it would seem desirable to transform the mode 
series to a form where the first term corresponds to the model of a flat earth and the succeeding 
terms are corrections for curvature. This approach has been used by Pekeris [21] and more 
recently by Koo and Katzin [22] in their investigations of microwave duct propagation. From 
their work, it may be shown that 



PJ— cos e) 

sin viT 



iH<?\kS nP )-± (£f m 2) (kS nP )+terms in (p/a)*, (p/a)«, etc, (12) 



where v J rl/2 = kS n p and p = ad. Ho 2) and II^ arc 4 Hankel functions of the second kind of order 
zero and two, respectively. When the great circle distance p is reasonably small compared to 
the earth's radius a, only the first term in the expansion need be retained. 

For convenience in what follows, it is convenient to express the field components as a 
ratio to the quantity 

E =i(r 1 /\)I(Is(e- ik o)!p (77^12()7r); 

E is the radiation field of the source at a distance p on a perfectly conducting ground. Thus, 
for both the source and the observer near the ground, it is not difficult to show^ that 



where 

where 

and 
where 



E 2 =WE a 






E,= -SE 



S^ 



T p 



e tt >T,5nS n m»(kS n p) 

?? = () 



H 4 ,= — TE () /r] 



T ^ - T !L e ^ S 5 n S n H{*\kS nP )= -N S S 



(13) 



(14) 



(15) 



where terms containing (p/a) 2 , (p/a)*, etc., have been neglected. 

When &p>>l, corresponding to the "far-zone," the above expressions may be simplified 
since the Hankel functions may be replaced by the first term of their asymptotic expansion. 
This leads to the compact result 



W 

s 

T 



Xp/x)%T 



(A/X) 



2ttp 




X 




71 = 



si 

-Sl/N s 

si _ 



c -*2,rS„p/\ 



(16) 



which is valid for p>>\. As expected, the ratio of W to T for a given mode is S n which for 
low order or grazing modes is of the order of unity. The ratio of S to T, quite generally, is 
— l/Nff which is very small compared to unity; in fact, it vanishes for a perfectly conducting 
ground as it must. 
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The excitation by a horizontal electric dipole may be treated in a similar manner. A com- 
plete discussion of the theory is given elsewhere [17]. The expression for the vertical field E z 
of the horizontal dipole Ids which is located at the origin and oriented in direction $ = may be 
written in the form 

E Z =SE cos $ (17) 

where S and E are defined above. By using the reciprocity theorem, this result can also be 
deduced directly from eq (15). Comparing the above formula for E z (for a horizontal dipole) 
with corresponding formula for E z , for a vertical dipole, it may be easily shown that 

E z (for horizontal dipole) cos0/ 1 \i cos 



E z (for verticle dipole) N g 



O-jW'-tc* (18 » 



being valid in the asymptotic or far zone. In terms of the ground conductivity <r g and dielectric 
constante^, 

1 __/ ie^oo \i^/iei(S\i 

Ng \<T g +ie g wj \ (Jg ) 

This is a small quantity. 

The above formulas for W, S, and Tare valid only if p< <a . If the first curvature correction 
term is included, it is a simple matter to show, in the "far zone", that this amounts to multi- 
plying the right-hand side by the factor 

Now, returning to the mode series for a spherical earth and using the far field approximation 
for the Legendre function, it turns out that IF, S, Thave precisely the same form as in eq (16) 
except the factor (0/sin 6)* occurs in the right-hand side. Noting that 

(f) \h ff2 

- — - ) ^1+-^+ terms in 4 , 6 , etc. 
sin 6/ 12 

it is apparent to a second order that the derived curvature correction is the same as from 
the Hankel function series. 

5. Distance and Frequency Dependence 

Certain features of the mode series are best illustrated by a graphical plot of field strength 
versus distance under various values of the parameters. Since N g ^>y>N t it may be safely 
assumed for ELF that A=l/7V f . Then on the assumption that the ionosphere is behaving 
like a conductor, it is possible to write 



Ni= 
where 



-0-f) ! (19, 

o)l _ (plasma frequency) 2 



collision frequency 



The above relation for N t is an approximation which is valid when the collision frequency 

v between electrons and neutral ions is much less than the angular frequency o>. Since v is of 

the order of 10 7 in the lower E region, this is certainly valid at ELF. The (angular) plasma 

frequency o> is given bv 

2 <?N 
w6= 

where iVis the electron density, e and m are charge and mass of the electron, and e is the dielectric 
constant of free space. co r , which is proportional to the ratio N/v, is believed to be of the order 
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of 10 5 in the lower region of the Slayer. At least this appears to be the effective value deduced 
from vlf observations for highly oblique incidence [23]. At ELF it appears that the waves 
are reflected at higher levels in the ionosphere and the effective value of co r is somewhat larger, 
being possibly of the order of 10 6 . 

It should be noted in passing that co r may be replaced by ajeo where cr f is the effective con- 
ductivity of the ionosphere and € is the dielectric constant of free space. 

To indicate the variation of the field as a function of distance, the magnitude and phase of 
W^and Tare shown in figures 1 to 4 for the typical values: co r = 5X10 5 and h = 90 km. Actually 
|W T | and \T\, both divided by the square root of the distance p (in kilometers), are plotted 
since the curves become linear at larger distances. The slope of these linear parts of the curves 
are proportional to the attenuation rate (in decibels per 1,000 km) as usually defined. It is 
seen immediately that at the shorter distances, the curves are no longer straight. This im- 
mediately indicates that considerable caution should be exercised in computing attenuation 
rates from spectral analyses of "sferics." For example, it is only when the distance p exceeds 
about one-sixth of a wavelength, is it permissible to assume that log (|M^|/\/p) or log (\E Z \ X\/p) 
vary in a linear mariner with distance p. Similarly, the phase of W (or E z ) varies in a linear 
manner only when the distance exceeds about one-half wavelength. Similar remarks apply to 
the magnetic field variations with distance. 




1000 l?00 1400 1600 1800 2000 

DISTANCE, p (km) 

Figure 1. The normalized magnitude of the electric 
field as a function of distance from the source for 
h = 90 km. 

(p is expressed in kilometers in both ordinate and abscissa). 




1000 1200 1400 1600 1800 2000 
DISTANCE, p(km) 

Figure 2. The normalized phase (lag), in radians, of the 
electric field as a function of distance for h = 90 km. 




800 1000 1200 1400 

DISTANCE, p (km) 

Figure 3. The normalized magnitude of the magnetic 
field as a function of distance from the source for h = 90 km. 




200 400 600 800 1000 1200 1400 1600 1800 

DISTANCE, p (km) 

Figure 4. The normalized phase (lag), in radians, of the 
magnetic field as a function of distance for h = 90 km. 
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At the extremely short distances where p is somewhat less than the ionospheric reflecting 
height (i.e., 90 km), the magnitude of W varies essentially as the inverse square of the distance 
and the phase approaches 180°. In this same region the quantity T varies as the inverse of 
the distance and the phase approaches 90°. The general characteristics of the field at these 
extremely short distances are in accord with the near field behavior of a dipole on a flat ground 
plane when the presence of the ionosphere is neglected. In this rather trivial case 



11 



and 



V h kV/ 



(20) 
(21) 



To shed further light on the nature of the ELF fields in the waveguide, the radial impedance 
of the wave is now considered. Noting that 



ET\ = W 

2 = 



(22) 



which is by definition the (normalized) impedance of the wave looking in the radial or p 
direction. First it should be noted that for &p>>l, 



and for kp << 1 and p << h, 



W 1 

— ~ 1 

T- 

W 1_ 

T~ikp 



For the intermediate and interesting range, kp is comparable to unity and p is somewhat greater 
than h. Using the numerical values of W and T mentioned above, the ratio IW/^I and the 
phase (lag) defined by arg T- arg W are plotted in figures 5 and 6 and as a function of distance 
for various frequencies. As before, co r =5X10 5 and h=90 km. 



lit 




J 1 I Mi l l 



100 200 300 500 700 1.000 

DISTANCE, p(krr\) 



Figure 5. The magnitude of the radial wave impedance. 




200 300 500 700 1000 2000 

DISTANCE, />(km) 

Figure 6. The phase (lag) of the radial\wave impedance. 
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The impedance ratio W/T is independent of the frequency spectrum of the source provided, 
of course, it may be represented by an equivalent vertical electric dipole. In fact, this complex 
ratio could be easily calculated from the frequency spectra of the waveforms of the vertical 
electric field and the horizontal magnetic field of an atmospheric. The observed variation of 
magnitude or phase of W/T as a function of frequency should then provide a basis for distance 
measuring. Such a scheme, while admittedly crude, requires only one receiving station 
equipped with a vertical whip and a loop antenna. 

6. Effect of the Earth's Magnetic Field 

In the preceding it has been tacitly assumed that the ionosphere is behaving as an isotropic 
homogeneous conductor. In this case, the refractive index N t of the ionosphere may be written 
[24] 

where 

_ oil _ (plasma frequency) 2 
7 v (collision frequency) 

When the earth's magnetic field is steeply dipping, it is appropriate to invoke the quasi-longi- 
tudinal approximation of Booker [24]. In this case the refractive index is double valued, one 
corresponds to the ordinary and the other the extraordinary, thus 



where 
and 



A^fl-i ^ exp (±ir)l * (24) 



_cu L _longitudinal component of gyrofrequency 
v collision frequency 



%^ 



(v 2 +u\)\ 



Using this model for the ionosphere, the reflection coefficient for a sharp boundary has been 
derived by Budden [25]. Adapting this result to ELF it has been shown [17] that the formula 
for S n has the same form as the isotropic medium if A is defined by 

A ^+(t) ' cos T/2 - (25) 

Since ti r =u r cos r it is possible to write the results in the same form as the isotropic ionosphere 
if we set 



where 



(m^^t] 



(Ni) ett and (w r ) eff are the effective values of N t and w r , respectively. 
Specifically, 



/ \ C0S T /o*\ 

(w ' )e »= W ' (^2?' (26) 

395 



In other words, the effective conductivity of the ionosphere is modified by the factor 
cos r/ (cos r/2) 2 which varies from unity to zero as r varies from to x/2. Thus, the attenuation 
is increased as a result of a steeply dipping (or vertical) magnetic field. In fact 



attenuation with magnetic field ^cos (r/2) 
attenuation without magnetic field = / cog \i 



(27) 



If r<0 (i.e., ul^v), this ratio becomes unity and the influence of the earth's magnetic 
field vanishes. At the level in the ionosphere where ELF waves are reflected, it is expected 
that o) and v are comparable, both being of the order of 10 6 . Assuming that they were actually 
equal, r becomes 45° and the ratio of the attenuation rates is 1.05. Thus, it is only when 
co L is greater than v does the earth's magnetic field appreciably influence the attenuation. 
The ratio is plotted in figure 7 as a function of tan r or co L /v. This illustrates the situation 
clearly. 




Figure 7. The influence of earth magnetic field on 
attenuation. 

The abscissa is the ratio of the vertical or longitudinal component of 
the gyrofrequency, col, and the collision frequency v. 



5 10 20 

to. /V OR tan r 



The preceding results are subject to the validity of Booker's quasi-longitudinal approxima- 
tion of the Appleton-Hartree equation [24]. The validity of this approximation requires that 



4co 2 a>2 ^ 



IV & <*>/ 



(28) 



where co L and co T are the longitudinal and transverse components of the (angular) gy rofrequency. 
Clearly this condition is violated when the transverse component of the earth's magnetic 
field is large such as for propagation along the magnetic equator. This case, however, has 
been considered by Barber and Crombie [26] who derived explicit results for the reflection 
coefficient at a sharply bounded ionosphere with a purely transverse magnetic field. Adapting 
their results to ELF, it is not difficult to show that A should now be replaced by [17] 



i= w+(B 



x 



(29) 



where 
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2 / ,2 , ,2 

COrp CO 



Now at ELF, co<Cco r and thus X is very close to unity since co T and v are of the same order 
of magnitude. 3 Consequently, a transverse magnetic field appears to have a negligible effect 
on the attenuation and the phase for ELF. 4 

A reasonable conclusion from the above is that only the vertical component of the earth's 
magnetic field is effective in ELF propagation. Furthermore, the ionosphere is effectively 
an isotropic conductor even in the presence of the earth's magnetic field. 

7. Effect of an Inhomogeneous Ionosphere 

The electron density in the actual ionosphere usually increases with height in the E region. 
Furthermore, the collision frequency decreases with height. Thus the effective value of the 
refractive index cannot be assumed constant. An approach is to let the refractive index 
increase or decrease from some initial value N at height // in a monotonic fashion. 5 Choosing 
an exponential variation the refractive index as a function of height is explicitly 

N(z)=hO forO<2<A 

= iVexp \{z-h)ll\ for z>h CM) ) 

where / is a scale factor. If />(), the refractive index is increasing with height and, if I 0, the 
refractive index is decreasing with height. It has been shown elsewhere [16] that the factor (J 
has t he forms 

q= MMe> iil>{) (31) 

K x (iNkl) 

and 

Q=s U-iNkl) ifz (32) 

Ixi—iNkl) 

where Io,Ii, K , and K x are modified Bessel functions with their conventional meaning. 

The preceding results are valid for |iV|>>l which is well satisfied at ELF. To the same 
approximation 

Nm[l-i^~f^y e-**» (33) 

where aJ r =coo/T> in terms of the plasma frequency oo and collision frequency v. It then follows 
that the imaginary and real parts of the propagation constant kS are given by 

and 

KekS =k+~ h (j^yP'(x) (35) 

where 

P(x) = |Q|sing-arga)V2 



P'(z)=!<2|co8 0-arg<2)V2 



3 KTote the positive values of cor correspond to propagation from east-to-west while negative values cot respond to propagation from west-to-east. 

4 If cot is somewhat greater than v, it is seen i hat to a fust order 

Xezi-ioT/p) (ico/co,)V2 
which will modify the Re -So to a slight extent but not Im So. 

5 In this section the influence of the earth's magnetic field is neglected. 



544328—60- 



397 



and 



c=\klN\ = \kl\(u r /a)* 



The functions P(x) and P'(x) approach unity if a* is sufficiently large. This limiting case cor- 
responds to the homogeneous ionosphere. Thus P(x) is the modification of the attenuation and 
P'(x) is the modification of the phase resulting from a nonhomogeneous ionosphere. It may be 
observed that the attenuation is generally lower if the refractive index is an increasing function 
with height. On the other hand, the attenuation is increased where the refractive index is 
decreasing with height. 

To illustrate the behavior of the attenuation rates as a function of frequency, they are 
plotted in figures 8a and 8b in terms of decibels per 1,000 km of path length for a frequency 
scale from 50 cps to 1.5 kc. Values of the scale distance I are indicated in kilometers. The 
height h of the lower edge of the ionosphere is taken to be 90 km. The values chosen for co r 
are 10 5 and 10° which are not inconsistent with previous work. 




Figure 8a. 



pr 



008 0.1 0.2 0.3 

FREQUENCY, kc 

Attenuation at ELF for an exponential 
ofile, \\~90 km and aJ r =10 r '. 




0.050.06 0.08 0.1 0.2 0.3 

FREQUENCY, kc 

CALCULATED 

MEASURED (CHAPMAN AND MACARIO) 

Figure 8b. Attenuation at ELF for an exponential 
profile, h = 90 km, Z T =10 b . 

The dashed curve corresponds to nighttime experimental data of 
Chapman and Macario. 



It is seen that the curve for l=co corresponding to a homogeneous ionosphere has a slope 
of 1/2 as it should. Depending on the sign of I and its magnitude, the slope may be modified 
considerably for the inhomogeneous ionosphere. The dotted curve in figure 8b corresponds to 
nighttime experimental data from Chapman and Macario [27] who obtained it from the spectral 
analyses of large numbers of atmospherics recorded in London. Clearly this experimental 
curve appears to fit fairly well on the l = 30-km curve, at least in the range 100 cps to 1.5 kc. 
There is no reason to expect any better fit than this because of the idealized profile assumed. 
Furthermore, the atmospherics analyzed by Chapman and Macario were often quite near the 
receiving location in terms of a wavelength. 

8. Propagation of ELF Pulses 

In the foregoing discussion, it has been assumed that the source is time harmonic. In most 
cases of practical interest, the current in the source dipole is of a transient nature such as a 
surge. The radiated fields are also transient in nature. While the resulting waveforms may 
be transformed to the frequency plane via spectral analyses, it is often more convenient to 
study the waveforms themselves [28 to 31]. 
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The source is assumed again to be equivalent to a vertical electric dipole on the earth's 
surface, but now its current-moment is a function of time and is denoted p(t) which is zero 
for£<C0. At a distance p along the earth's surface, the resulting vertical electric field is e x {t) 
and the horizontal magnetic field is h^{t). The Laplace Transform of the current moment is 
denoted Pq(s) or more explicitly 



P (8) = Lp(t)=\~ e-'p(t)dt 



(36) 



where s is the transform variable and may be formally identified 
with io). Also 

E g (s)=Le z (t) 
and 

where L is the Laplace-Transform operator defined above. Now it is assumed thai the im- 
portant frequencies in the spectra are sufficiently low that only the zero-order mode need be 
retained. Thus 



where 



s '=N(£)j" 



(37) 



and where K and K x are modified Bessel functions of argument sS p/c. These may be asymp- 
totically represented by 



and 

for large x, where 
Thus 

E.(«) 



^KSH i+ t+§°+ • • •] 

A = -1/8, £o=9/128, ^ = 3/8, B x =- 15/128. 



(38) 



(39) 



E.(«) „j^(2\h(c_\i i 



i £,ft/cY 



■^G) 2 " 1 



8 p S 2 \pJ 



(40) 



where x has been replaced by sS p/c in the exponent but is replaced by sp/c elsewhere. This is 
justified since S is near unity. In fact, the binomial expansion 
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may be truncated beyond the second term in most cases of interest, although the third term is 
retained since it doesn't complicate matters. 

The current moment is now taken to be an impulse, that is 

p(t)=p»8(t) 
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where 8(t) is the unit impulse function. Thus 
The transforms may now be readily inverted 

(42) 
where 

P a (7 , )=2(iJ /2 exp(- i L) (44) 

P 6 (7 , )=4(IJ /2 exp (~^)-2 ^ erfc (^) (45) 

t>=t- P ,c and *-£(*)* 



and 



It is immediately seen that if (ct'/p)«lj the field responses are both proportional to P(T). 
To illustrate the nature of this function, it is plotted in figure 9 as a function of T. To facili- 
tate the application a multiple scale is included which relates actual time V in microseconds 
with the parameter T f or distances p of 1,000 to 4,000 km. The height of the ionospheric 
reflecting layer is taken as 70 and 90 km. The quantity B indicated on the curves is related to 
the ionospheric conductivity by 

r>_ 1 __2tc 
60cjA <j) T h 

The values of B shown, namely 0.05 and 0.1, are typical. For A=90 km, these correspond to 
w r values of 4X10 5 and 2X10 5 , respectively. 

The responses e z {t) and h^t) to a general source p(t) can be expressed in terms of the im- 
pulse responses e° z (t) and hl(t), respectively, by using the convolution theorem. Thus 

^W=4 f P(t-T)e%T)dr-u(t-p/c) (46) 

Po J P/c 

and similarly for h^if). For example, if a special analytical form for p(t) is chosen such as 

*>=f(^) ! >K-t> <«> 

the convolution integral may be evaluated to give 

«t)*-*®*J^gf fa «p[g(;)]P(J'/r) (48) 

where 

The function P(T) was defined above and is the same as the one plotted in figure 9. Now, 
however, T is to be identified with t'/y rather than t'/p. The special form for p(t) given above 
is a unidirectional pulse which reaches its maximum value at t=t m and decays to 7 percent of 
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Figure 9. Transient response of the zero-order mode for 
an impulsive vertical dipole source. 

The lower chart facilitates conversion from the parameter T to actual 
timet'. 

its maximum value at t=St m . If the duration, U my of the source pulse is much less than 
the response e z {t) approaches e\t). This is not surprising since p\t) approaches p 8(t) if t m 
approaches zero. 

The general effect of the finite duration of the source pulse is seen to "stretch" the waveform 
of the electric and magnetic field. This effect has been discussed at some length in a recent paper 
by the author entitled "On the theory of the slow-tail portion of atmospheric waveforms." This 
is to appear in the July issue of the Journal of Geophysical Research. 

While most cloud-to-ground lightning strokes may be represented by a vertical electric 
dipole, it is believed that cloud-to-cloud discharges are better represented by a horizontal 
electric dipole [32]. In the ELF portion of the spectrum the height of the discharge is very 
small compared with the wavelength, thus eq (17) is valid. If the source moment is again 
p(t) then the transform for the vertical electric field, in the far zone, is given by 



^m® 



sPMe-'W 



w 



COS 



where </> is the angle subtended by the horizontal dipole and the direction to the observer. 
Letting the source be an impulse [i.e., p(t)=p 8(t)] } the transform may be readily inverted to 
give 



M"&(£ffk($<* T > 



(49) 
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where 



and where, as before, 



Q(T)=~ (j^-3) exp (-^} (50) 



T=t>,fi and t-frffi 

An immediate generalization is a source which may be imagined as an inclined electric dipole. 
Lightning discharges from cloud-to-ground and cloud-to-cloud would be seldom purely vertical 
or horizontal. In view of the small dimensions of the discharge paths in terms of a wave- 
length, it is permissible to replace the inclined channel by superimposed vertical and hori- 
zontal electric dipoles. The resulting responses of the radiated field are thus obtained by super- 
position. For example, if the current moment of the vertical electric dipole is p 8(t) then the 
current moment of the horizontal electric dipole is gpo5(t) where g is a positive or negative 
dimensionless number. For lightning strokes with long horizontal sections, g could be large 
compared to unity [32]. As before, the direction 0=0 corresponds to the direction of the 
horizontal dipole component. 

Invoking the above mentioned principle of superposition, it readily follows that the 
far zone expression for the inclined dipole is given by 



.(J)sg(0 [P(T)+GQ(TW-^ 



where 



\<W P 



COS <f>. 



To indicate the bahavior of the waveform, the quantity 

P(T)+GQ(T) 



S(T)- 



1+101 



is plotted in figure 10a and 10b for both positive and negative values of G. The quantity 
S(T) characterizes the waveform of the inclined dipole energized by an impulsive current. 
The special case 6?=0, such that S(T) = P(T), corresponds to a vertical dipole source, whereas 
the special case G= 00 corresponds to a purely horizontal dipole source. Such a range of G 
values can be expected since cloud-to-ground and intra-cloud discharges are primarily vertical, 
whereas the cloud-to-cloud discharges and air flashes are mainly horizontal and may have 
either polarity. 
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Figure 10. Transient response of the zero-order mode for an inclined dipole source. 
When G=0 the dipole is vertical and when G=±<*> the dipole is horizontal. 
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The modification of the pulse shape of the ELF waveforms, as a result of the inclination of 
the current channel, would appear to be an important factor in interpreting observed data. 
For example, Pierce [9], Lieberman [29], Tepley [31], and Jean [33] have all observed ELF 
waveforms which are strikingly similar to those shown in figures 10a and 10b. For example, 
the observed reception of a pulse which has a second half-cycle of the same magnitude as the 
first half-cycle can only be reconciled with an inclined source with G equal to about ±0.2. For 
certain values of G (e.g., small negative values), a third half-cycle of relatively small amplitude 
is also produced. 

9. Concluding Remarks 

In the present theoretical study it lias been assumed that the ELF signal propagates entirety 
in the space between the earth and a concentric ionosphere. It is quite possible that whistler 
type propagation may also be important. Furthermore, sources of ELF energy may be present 
in the exosphere due to the possible existence of ionized hydrogen whose gyrofrequency is of the 
order of 900 cps [34]. Other types of ions such as sodium may also be significant as pointed out 
by Aarons [14]. Despite these extraneous factors, the evidence that the "slow tails" are propa- 
gated in the earth-ionosphere waveguide is overwhelming. 

It is the opinion of the author that the only really questionable assumption is the neglect of 
heavy ions on the constitutive properties of the lower edge of the E region or the top of the D 
region. Normally, it would be expected that the constitutive properties of the ionosphere 
would only be influenced by the electrons because of the extremely large ratio between the 
masses (i.e., of the order of 10 7 ). However, if the ratio of the number of heavy ions to the 
number of electrons is of the order of 10 3 or 10 4 , as has been suggested by Scott [35] for the 
Arctic Slayer, flic influence of heavy ions may be significant at frequencies less than 100 cps. 
In particular, it should be noted that the gyrofrequency of the sodium ion is approximately 30 
cps. However, if the operating frequency is somewhat greater than the gyrofrequency of heavy 
ions, their influence is expected to be minor. 



The author is indebted to Mrs. Nancy Carter and Mrs. Alyce Conda who carried out many 
of the calculations, to Mrs. Patricia Murdock who typed the manuscript, and to Mrs. Barbara 
Bolton who prepared the illustrations. 
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